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Abstract 

We investigate the tensor and the scalar perturbations in the symmetric bouncing universe 
driven by one ordinary field and its Lee- Wick partner field which is a ghost. We obtain 
the even- and the odd-mode functions of the tensor perturbation in the matter-dominated 
regime. The tensor perturbation grows in time during the contracting phase of the Universe, 
and decays during the expanding phase. The power spectrum for the tensor perturbation 
is evaluated and the spectral index is given by nx = 6. We add the analysis on the scalar 
perturbation by inspecting the even- and the odd-mode functions in the matter-dominated 
regime, which was studied numerically in our previous work. We conclude that the comoving 
curvature by the scalar perturbation is constant in the super-horizon scale and starts to 
decay in the far sub-horizon scale while the Universe expands. 
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1 Introduction 

Inspired by the original work [1] of Lee and Wick on the quantum electrodynamics with higher- 
derivative propagators, Grinstein et al proposed the Lee- Wick standard model [2] . In this model, 
the minimal set of the higher-derivative field was introduced, and the corresponding N = 2 Lee- 
Wick (LW) formalism was obtained. The N = 2 LW theory was extended to the = 3 one in 
Refs. |3]. For the higher-derivative model consisted of one scalar field, the general LW formalism 
was constructed in Ref. [1]. (See also Refs. [SHZ] for recent works for the LW field theory.) 

The N = 2 scalar-field LW model consists of one ordinary scalar field and its LW partner field 
which is a ghost. When this model is applied to cosmology, the background universe experiences 
bouncing due to the energy-condition violating ghost field. For the bouncing-universe models, for 
example the ekpyrotic model [8], the perturbation becomes singular in general about the bouncing 
point. However, very recently the N = 2 scalar-field LW model was investigated in Ref. [9], and 
the authors found that the density perturbation at the bouncing point is nonsingular. (See also 
Refs. [ID]-) Although the stability of the bouncing universe about the anisotropic perturbation is 
still questioned, at the linear level it is believed that the LW bouncing universe is stable about 
the perturbation. (See also Refs. [TT] for other types of nonsingular bouncing universe models.) 

In Ref. [12], we investigated the scalar perturbation of the symmetric Lee- Wick bouncing 
universe, in which a new type of initial vacuum solution was discovered. We discussed the growth 
of the initial perturbation and the resulting late-time power spectrum by introducing the even- and 
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odd- mode analysis. In this paper, we investigate the tensor perturbation applying the even/odd- 
mode analysis, and complete the discussion on the scalar perturbation which was studied in our 
previous work. 

The paper is organized as following. In Sec. 2, we review the set-up of the N = 2 symmetric 
Lee- Wick model and the evolution of the background universe. In Sec. 3, we solve the tensor 
perturbation, and obtain the power spectrum and the spectral index. In Sec. 4, we extend the 
mode analysis to the scalar perturbation, and comment on the shortness of the analysis. In Sec. 5, 
we conclude. 



2 Symmetric Bouncing Universe 



The N = 2 Lee- Wick model consists of one ordinary scalar field ipi and one ghost field ip2- The 
Einstein-Hilbert action with these Lee- Wick matter fields is given by 

R 



LW 



16nG 



n=l 



(2.1) 



where we assumed no interacting potential between the matter fields. For > mf, the ghost 
field plays an important role about the bouncing point, and becomes subdominant for the rest of 
time. The metric ansatz for the isotropic fiat Friedmann- Robert son- Walker universe is given by 



ds"^ = -dt^ + a^{t)dx'dx\ 
The Friedmann equations and the scalar-field equations are then given by 

n=l 
2 

H = -47rG^( 



f2.2) 



n=l 



3if 0„ + m^(y9„ = 0, (ra = l,2). 



(2.3) 

(2.4) 
(2.5) 



The ghost field violates the null energy condition, and thus the background evolution can have 
a bouncing from contracting to expanding. 

When the Universe undergoes bouncing (let us set the bouncing moment at t = 0), the velocity 
of the scale factor vanishes, d(0) = 0, i.e., the Hubble parameter becomes -ff(O) = 0. From the 
Friedmann equation (12. 3p . we have 



a(0) = ^ [Pi-^l + ^Wl 



2 21 



t=0 



0. 



(2.6) 
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For the symmetric bouncing, the scalar field need be symmetric (even or odd) about t = as 
well as the scale factor. In order for to be symmetric, the necessary condition is 



There are four cases satisfying the above symmetric condition for as it was studied in Ref. |12j . 
Among them, we selected most conventional conditions for the bouncing universe. 



which indicates that <y9„'s are even functions. This conditions are depicted in Fig. 1. At t = 
the ordinary field rolls down the potential from rest, while the ghost field climbs up the potential. 
Afterwards, the two fields oscillate about = behaving as dusts. Then the background universe 
undergoes the matter-dominated expansion. Due to the background expansion, the matter-field 
oscillations are damped. Setting the ghost field more massive, m\ < m^, the late-time evolution 
is dominated by the light ordinary field. The background fields at late times approximate as 



where a and 0i are a numerical constant. Since the background fields a{t) and (fn{t) are symmetric 
about t = 0, starting from some initial moment t = ti < 0, the scalar fields oscillate about ipn = 0, 
and the scale factor decreases, i.e., the Universe contracts. Meanwhile, the amplitude of (pn 
increases and reaches the maximum at t = 0, and then evolves as stated above for t > 0. (The 
numerical plot of the background fields are plotted in Fig. 1 in Ref. |12j.) 

We consider linear perturbations in this background. The initial scalar and tensor perturba- 
tions are produced during the contracting phase at t = ^ 0. The initial perturbation for a 
given wave number crosses the horizon four times in total. The horizon scale \H~^\ is plotted in 
Fig. 2. The initial perturbations are produced deep inside the horizon at t = ti, and then cross 
out the horizon. When the perturbations finally cross in the horizon during the expanding phase 
at t > 0, the power spectrum is evaluated. Afterwards, they enter the nonlinear regime, which is 
beyond our concern in this work. 

For the scalar perturbation in the symmetric background universe, rather than solving the 
perturbation from t = numerically, one is allowed to solve from t = to either direction owing 
to the symmetry. In Ref. [12], we analyzed the perturbation about t = 0, and found that the 
perturbation is composed of two linearly independent modes (even and odd modes). Then we 
solved numerically the even and the odd mode from t = in the t > region, which is numerically 
stable. Due to symmetry, we need not solve the perturbation for t < 0. This is the advantage of 
solving the even and the odd mode rather than other modes in numerical calculations. 



</.„(0) = 0, or V9„(0) = 0. 



(2.7) 



0^(0) = and m?(^?(0) = mlipl{Q) ^ 0, 





(2.9) 



4 



3 Tensor Perturbation 



In this section, we investigate the tensor perturbation in our symmetric bouncing background. The 
tensor perturbation is decoupled from the scalar and the vector ones in the linear perturbation 
theory. In order to describe the tensor perturbation, let us consider the perturbed metric in the 
conformal time. 



ds^ = a^{r]) [—dr]'^ + {5ij + hij)dx'^dx-'~\ , 
where the tensor perturbation hij satisfies the gauge conditions. 



dihij = 0, ha = 0. 



(3.10) 



(3.11) 



With these conditions, the tensor mode hij possesses two degrees of freedom corresponding to two 
polarizations of gravitational waves, and can be expanded as 



A=+ - 



(3.12) 



where efj represents the polarization tensor. Inserting ( I3.10p and (13.12^ into the Einstein-Hilbert 
action, we obtain the action for the tensor perturbation, 



A=+, 



S't = - E / dr]d^x [drifix] ^ - [difix) 



2 a _2' 
a 



(3.13) 



where ' = d/drj and we used the relation e^je'^, = 6y. The general solution to the equation of 
motion for fi\ is given by 



d^k 
(27r)3/2 



Ak-x 



(3.14) 



where Hxiji] k) satisfies 



/"A + - - = 



(3.15) 



with the normalization for a given polarization mode A, 



(3.16) 
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3.1 Initial and Super-Horizon Scale Perturbations 

Let us investigate the mode of the tensor perturbation, which is produced in the matter-dominated 
regime during the contracting phase. The vacuum states of small comoving-wave numbers (long- 
wave lengths) exit the Hubble horizon while the Universe contracts. After they exit the Hubble 
horizon, they experience the bouncing phase. At the bouncing point, the Hubble radius becomes 
infinity as shown in Fig. 2. Therefore, all the Fourier modes in the super-horizon scale re-enter into 
the Hubble horizon near the bouncing point, and then re-exit the horizon shortly after bouncing 
during the expanding phase. 

When the background universe is in the matter-dominated regime, the scale factor becomes 
a ^ at^/^. One can transform the cosmological time to the conformal time by dt = adr], then 
obtain t = and a ^ a^rj'^/d. Then one has a"/a ~ 2/?]^, and the equation f l3.15p becomes 

/^A + - ^) -"A ^ 0. (3.17) 

This equation looks the same with that of the single-field inflationary models, but the perturbation 
described by the tensor mode hij is different because the scale fact a is different. 
The general solution to Eq. fl3.17p is given by 



f^xiv,k) = \B^{k)H^^\-kr^) + Ax{k)Hf\-kr,)\ , (3.18) 

V Z L 2 2 . 

= + (3.19) 

= l^f(ri) + (3.20) 

where Hn''^^ represents the Hankel functions, and \B\\'^ — \Ax\^ = l/2k from the normalization 
fl3.16p . If we introduce the even and the odd mode which are linearly independent real functions, 

p ,, , sin(A;?7) ^ . . cosikri) 

fif = cos(/c77) //O = sm{kv) + — (3.21) 

fi^^ and /i^™ are expressed by a linear combination of these even and odd modes. By setting 
Ax = Oi + B\ = Pi + i(32, one obtains the following relation. 




(3.22) 



(i) Initial Perturbation: The initial perturbation is produced deep inside the horizon when the 
Universe contracts. Then the perturbation crosses out the horizon when its physical wave-length 
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scale is comparable with the horizon scale, 

Xp,, = ^^\H-'\ =^ \kr]\r^2. (3.23) 
In the far sub-horizon limit, \kT]\ ^ 2, the even and the odd modes are approximated by 

fif ^ cos(A;?7), fix ~ sin(A;?7). (3.24) 
We adopt the Bunch-Davies vacuum for the production of the initial perturbation by taking only 



the positive-energy mode. Ax = (a^ = ^2 
Eq. (I3.19p . we obtain the initial perturbation. 



0). Then using Eq. (]3.24p . or directly from 



2k 



(3.25) 



where 9 = tan'^i^^ / (3^) is a phase factor from the complexity of Bx ■ The amplitude of the tensor 
mode becomes 



IT I If^M 1 
\hij \ oc oc — , 

a 7]-^ 



(3.26) 



so it grows initially in the contracting universe. 



(ii) Super-Horizon Scale Perturbation and Power Spectrum: Let us discuss the pertur- 
bation in the super-horizon scale during the expanding phase, and its horizon crossing. We have 
|A;?7| <^ 2 in the super-horizon scale, and I/ct^I ~ 2 at the horizon crossing. Therefore, the second 
term in in Eq. fl3.2ip is dominant, or comparable to the first term. The solutions yuf and 
in Eq. f l3.2ip are valid once the background is matter-dominated. Therefore, the conditions 
imposed for the initial perturbation such as the Bunch-Davies vacuum and normalization, are still 
valid in the super-horizon scales. 



A, 



0, 



and \Bx\' = {P^^f + {P^) 



2k' 



(3.27) 



Then from the solutions /if and //^ in Eq. f l3.2ip . we have the relation. 



/^A/iA 



.E\2 



.0\2 



1 

2k 



1 + 



{krif 



(3.28) 



Using this result we can obtain the power spectrum. 



Pi 



/^a| 



27r2 



1296GP 



1 + 



{krjy 



(3.29) 
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The power spectrum of the tensor perturbation decays as the Universe expandsjl] The power 
spectrum at the horizon-crossing (|fc?7| ~ 2) is given by 



and the tensor spectral index becomes, 



3.2 Perturbation about the Bouncing Point 

Let us discuss the tensor perturbation about the bouncing point. The background fields, ^Pnif), 
a{t), and H{t) had been obtained in series expansion about t = in Ref. [12]. Using the scale 
factor 

a(t) = 1 + ]ht^ + lhr,f + ■■■ , (3.32) 
4 

where /i^'s are constants determined in Ref. [I2], we can transform the time coordinate t to rj 
about the bouncing point, 

t = v + ^v' + ---- (3.33) 



Then we have the scale factor in conformal time 

i^^^^ ■ 6 



a = 1 + ^ W + IKv"" ■■■ , (3.34) 



and also have 

a" 

— = 3/13//' + 5hr]^ + ■■■ . (3.35) 

a 

With this, the field equation fl3.15p for the tensor perturbation can be solved. Since a" /a is an 
even function of 77, the general solution can again be expressed by a linear combination of the even 
and the odd mode which are real functions, 



1/U 5U4(J 



(3.36) 



^The power spectrum Vt evaluated in Eq. ()3.29p is valid for the whole period of the matter-dominated regime, re- 
gardless of the background phase (contracting/expanding) and the perturbation scale (sub-horizon/super-horizon). 
[For example, when \kii\ ^ 2, it reduces to that of the initial perturbation approximated by Eq. p.26p .] Therefore, 
apart from the bouncing point, the amplitude of the tensor perturbation grows in time during the contracting 
phase, and damps in the expanding phase. 
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where 60 and bi are integration constants. 

In the previous subsection, we obtained the even and odd modes in the matter-dominated 
regime. For a given Fourier mode, the mode coefficients are constrained by the conditions imposed 
at the initial state. Once the coefficients are fixed, they do not change during evolution since 
they are functions of the Fourier mode k only. One can see easily that the exactly same linear 
combination of /if and /i^ in Eq. fl3.36p as in Eq. fl3.22p used for the initial perturbation, satisfies 
the normalization (I3.16P with bo = bi = \/2k in the limit of 77 — )■ 0. Therefore, one can say that 
the even (odd) mode in Eq. (13.360 obtained near the bouncing point, is continuously connected to 
the even (odd) mode in Eq. (I3.2ip obtained in the matter-dominated regime. In other words, the 
evenness (oddness) is preserved in the whole range. This fact will be useful when one is to solve 
the field equation (I3.15p numerically from the bouncing point using Eq. (13.361) . When one solves 
the field equation numerically from the bouncing point with the even (odd) boundary conditions 
at = using Eq. (13. 36 p . the numerical solution at large |?7| should reproduce the even (odd) 
solution in the matter-dominated regime. This story is valid also for the scalar perturbation in 
the next section. 

4 Scalar Perturbation 
4.1 Initial Perturbation 

In Ref. [12], we investigated the scalar perturbation. We solved field equations numerically, and 
analyzed the solutions using the even- /odd-mode technique. In this section, we complete the 
analysis of the scalar perturbation in a similar way with the tensor perturbation analyzed in the 
previous section. 

In order to discuss the scalar perturbation, we introduce the Sasaki-Mukhanov variable Q„ 
which reduces to the perturbation of the scalar field in the spatially flat gauge, Qn = Sipn- (For 
details of the scalar perturbation, see Ref. [I2]-) We introduce the field v = aQi for canonical 
quantization, where Qi = 6ipi corresponding to (pi. (For our model, the ordinary field ipi is 
dominant for the scalar perturbation, so one may consider v only.) The action is given by 



where z = a(p[/'H with "H = a' /a. The general solution of the equation of motion is given by 




(4.37) 



v{r], x) 



I 



v{r];k)af: + v*{T];k)a^_j. e 



ik-x 




(27r)3/2 
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where ^(77; A;) satisfies 



and is normalized as 



v" + k 



z" 

z 



w = 0, 



VV — V V = I. 



(4.39) 



(4.40) 



In the matter-dominated regime, the field equation fl4.39p can be approximated by 



V"{7]) + 



2 m?a^ 4 4mia^ / 2mia^ o 
-?7 - ' 



-7] sm 



27 



?7'' + 2a] 



v{r]) ^ 0. 



(4.41) 



In the far sub-horizon limit, |/c?7| ^ 2, the equation f l4.4ip can be further approximated by 



d{kriY 81A;6 
The solution to this equation is given by 



{kr]Yv{T]) ^ 0. 



v{ri) 



-TTT] 
~\2 



(1) 



rriia n 



(4.42) 



(4.43) 



where \A2\ 
v{ri) ^ 



y4ip = 1 from normalization (14. 40 p . For rj <^0, this solution is approximated by 

(4.44) 
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— < Ai exp 



— — V + - 
27 ' 3 



9 



- (Ai + A2 



cos 



rriia o 



-Ai + A2) 



sm 



mia o 
27 ' 



A2 exp 
— -\/3sin 
- -\/3cos 



27 ' 3 



rriia o 



27 ' 



If we introduce the even and the odd mode which are linearly independent real functions. 
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1 

- cos 



irriia-' rj 



9 



1 

— sin 



27 ^ 



(4.45) 
(4.46) 

(4.47) 

f and w^™ can be expressed by a linear combination of these even and odd modes. By setting 
^1 = /3i + ^/32 and A2 = ol\^ ia2, we obtain the following transformation, 

-a/3/3i + 132- V^ai - a2 - ^3/32 - «i + VSa2 

-f3i — V3(32 — «! + \/3a2 V^(3i — (32 — y/Sai — 0L2 

For the initial perturbation, we adopt the positive-energy mode, A2 = 0, i.e., ai = 0^2 = 0, then 
we have 







1 ylm 1 







f \ 


) 




e., ai = ( 



(4.48) 



,Rc 



Im 



-v^/3i + /32 -A-V3/32 
-f3i - 73/^2 V^(3i - (32 



,0 



(4.49) 
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4.2 Super-Horizon Scale Perturbation 



In Ref. [12], we observed from numerical calculations that the scalar perturbation gives rise to 
constant comoving curvature TZ in super-horizon scales. In solving the field equation fl4.4ip ana- 
lytically, it is not very possible to consider all terms in the equation. For this reason, it is difficult 
to analyze the super-horizon scale perturbations in terms of v because they are determined by 
the effect of mixture of terms in Eq. fl4.4ip . Instead, it is better to analyze the super-horizon 
perturbations in terms of the Sasaki-Mukhanov variable Qn{t)- The equation of motion for Qn is 
given by (see Ref. [12]) 



^2 SttG ^ 

Qn + 3ifg„ + — Q„ + mlQr, - ^ ^( 

1=1 



d f a 



dt\H 



(4.50) 



When the background is relaxed to the matter-dominate universe, we can approximate the equa- 
tion for Qi which is dominant as 



Q1 + -Q1 + 

V 



a 



2^4/3 



2 4mi 
+ sm(2mit + 20i^ 



(4.51) 



The constant comoving curvature at the super-horizon scale is achieved mainly when the m^-term 
is dominant. 

If we transform this solution to the canonical ffeld v = v{ri) in conformal time, we get 



v{r]) = aQi 



Ai exp 



mi a o 
27 ' 



. , .mitt'' o 
^2 exp ( ^^^V 



(Ai + A2) COS ( 
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27 



-A, + A,) sm ( -^T]^ 



(4.52) 



(4.53) 

(4.54) 
(4.55) 



If we introduce the even and the odd mode which are linearly independent 

.,3 



,p a^T]'^ I mia n 

V = — : — COS I — T^f] 



9 



27 



s 2 / 3 
o tt ?7 . / mia o 
V = — z — sm I 77 



9 



27 



(4.56) 



'0^'' and -u^™ are expressed by a linear combination of these even and odd modes. By setting 
^1 = Pi + and A2 = di + ia2, one obtains the following relation. 



/3i + di P2- «2 
P2 + ^2 + di 



(4.57) 
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In the previous subsection for the initial perturbation, we adopted the positive-energy mode. Once 
we pick up this mode, the positivity of the mode should remain unchanged even in the super- 
horizon scale. Therefore, we select the positive mode from the approximated solution v in the 
super- horizon scale by setting A2 = (di = 62 = 0). Then we finally get the transformation for 
the approximated super-horizon solutions. 



£,Im 




(4.58) 



The approximated solution w(?7) describes the perturbation at the super- horizon scale for 77 > 
as well as for 77 < 0. During the contracting phase, the initial perturbation is produced deep 
inside the horizon at 77 < 0, which is described by the solution f (77) in Eq. f l4.49p . Afterwards, the 
perturbation crosses out the horizon and becomes the super-horizon scale. Then, it crosses in the 
horizon and passes the bouncing moment. During the expanding phase, the perturbation crosses 
out the horizon soon after the bouncing moment and becomes super-horizon scale. Meanwhile the 
background universe relaxes to the matter-dominated universe, so the super-horizon perturbation 
is approximated by 0(77) in Eq. (14.581) . The power spectrum that we observe, is evaluated by this 
perturbation when it crosses back into the horizon. Afterwards, the perturbation starts to damp 
in principle, and is described by f (7/) in Eq. (I4.49P again deep inside the horizon. However, after 
the perturbation crosses the horizon, it reaches the nonlinear regime which 

The approximated super-horizon solution {"(t/) in Eq. (14.541) is not normalizable, i.e., it does 
not satisfy the normalization condition (I4.40p . The correct normalizable super- horizon solution 
v{ri) must require some corrections to the approximated solution. However, it is not very possible 
to obtain v{ri) in our situation. Therefore, let us evaluate the power spectrum in approximation 
using V, which will imply underlying physics to some extent. 

The power spectrum for the scalar perturbation is given by 

Vn = |^|7^|^ (4.59) 
where the comoving curvature is defined by 

n = i^iQi - ip2Q2) = fiQi - /2Q2 ~ /iQi. (4.60) 



^The perturbation entering deep inside the horizon at > is mathematicaUy described by the same solution 
v(ri) produced deep inside the horizon at rj < 0. Therefore, if one solves the field equation numerically for 
t > 0, one should get the solution transiting from v (which provides constant comoving curvature TZ) to v (which 
provides decaying TZ). In this sense, we misinterpreted in Ref. |12| where we commented that the shape of v 
(and corresponding T^^onst-j niaintains even after entering the horizon, while the behavior of v (and corresponding 
-^decaying ^ ^-i-j sub-dominantly implied in v in the sub-horizon scales. 
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In the last step, we used the fact that Qi is dominant. When the background is relaxed to the 
matter-dominated universe, using Eq. f l2.9p we have 



/i 



H 



The approximated comoving curvature is the given by 

IZ = — V, 



16nG 



■ cos 



(4.61) 



(4.62) 



then the approximated power spectrum which is averaged over time, becomes 



27r2 



27r2 a2 



Im\ 2 



16G, 
~37r 



(4.63) 



where B'^{k) = /^^ + /3| and the cosine factor has been averaged out. Unfortunately since v does 
not satisfy the normalization condition f l4.40p . we do not know the value of /3j. In order to provide 
the scale-invariant power spectrum, B'^{k) oc k~^ is requiredlfl 



5 Conclusions 

In this paper, we studied the tensor and scalar perturbations of the symmetric Lee- Wick bouncing 
universe. We adopted the even-/odd-mode analysis instead of decaying-/constant-mode analysis 
in the literature [9l[T3l[T3]. The even and the odd mode are linearly independent modes of the 
perturbation. The analytic solutions for these modes are obtained both in the bouncing point and 
the matter-dominated regime. The evenness and the oddness are preserved in these two regimes. 

For the tensor perturbation, the mode functions are analytically solved when the background 
is in the matter-dominated universe. These mode functions are valid both for the sub-horizon and 
the super-horizon scales. The initial tensor perturbation produced deep inside the horizon during 
the contracting phase at t < 0, grows in time. After the tensor mode crosses out the horizon 
during the contracting phase, it grows more rapidly. After the bouncing, the perturbation evolves 
in a symmetric way; it decays as the Universe expands. In the super-horizon scales during the 
expanding phase, the power spectrum Vt decays in time. Finally when it crosses the horizon, the 
spectral index is evaluated as nT = 6. 

■^If we consider \R.\'^ for far sub-horizon scales, /i is the same with Eq. (|4.6ip . and we can use the solutions vf^'^ 
and in Eq. (|4.49p evaluated by and zP in Eq. (j4.47p . Then we have 

^- (4-64) 
Therefore, {\'R\^) grows in the contracting phase, and decays in the expanding phase. 
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For the scalar perturbation, the initial perturbation obtained in Ref. [12] provides growing 
comoving curvature TZ. In the super-horizon scales, we obtained the approximated solutions 
for the even and odd modes. These approximated solutions provide constant comoving curvature. 
However, they are not complete enough to evaluate the fc-dependence of the power spectrum. (The 
approximated mode solution is not normalizable, so the mode coefficients which depend on k are 
not determined. The correct mode solution would require some corrections to this approximate 
solutions, but is difficult to obtain for our model.) Therefore, the scale invariance could not be 
examined. 

When the scalar perturbation finally crosses in the horizon during the expanding phase, the 
modes are of the same functional form with the initial perturbation due to symmetry. The cor- 
responding comoving curvature will decay, and the perturbation evolves finally into the nonlinear 
regime which is beyond our scope. As a whole, the comoving curvature by the scalar perturbation 
produced initially in the sub-horizon at t < 0, grows to become constant in the super-horizon. 
During the expanding phase at t > 0, the curvature remains constant in the super-horizon, and 
starts to decay after the perturbation crosses into the horizon. 
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Figure 1: At the bouncing point {t = 0), the ordinary field ipi and the ghost field ip2 are at rest, 
(pn = , with the values satisfying m\ip\ = rn^ip^. Afterwards at t > 0, two fields approach the 
center of the potential, and start to oscillate about ipn = 0. For t < 0, the behavior of ipn is evenly 
symmetric. For t > 0, the background universe expands due to two scalar fields, approaching the 
matter-dominated universe. For t < 0, the universe contracts in a symmetrical way. 
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Figure 2: Plot of the horizon scale |-f/^~ |. The perturbation for a given wave length is initially 
produced inside the horizon during the contracting phase at t < 0, and crosses the horizon four 
times. 
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